Devil's staircases and supersolids in a one-dimensional dipolar Bose gas 
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We consider a single-component gas of dipolar bosons confined in a one- dimensional optical lattice, 
where the dipoles are aligned such that the long-ranged dipolar interactions are maximally repulsive. 
In the limit of zero inter-site hopping and sufficiently large on-site interaction, the phase diagram is 
a complete devil's staircase for filling fractions between and 1, wherein every commensurate state 
at a rational filling is stable over a finite interval in chemical potential. We perturb away from this 
limit in two experimentally motivated directions involving the addition of hopping and a reduction 
of the onsite interaction. The addition of hopping alone yields a phase diagram, which we compute 
in perturbation theory in the hopping, where the commensurate Mott phases now compete with the 
superfluid. Further softening of the onsite interaction yields alternative commensurate states with 
double occupancies which can form a staircase of their own, as well as one-dimensional "supersolids" 
which simultaneously exhibit discrete broken symmetries and superfiuidity. 



The unprecedented control over experimental parame- 
ters in trapped ultracold atomic gases has substantially 
widened the range of possible exotic phases of matter that 
can be explored. Optical lattices can be used to simulate 
simple lattice models and/or vary the system dimension- 
ality, while the interatomic interactions can be varied via 
a magnetically-tunable Fcshbach resonance Thus far, 
the focus has largely been on contact interactions, since 
these generally provide a good description of atom-atom 
scattering in the low energy limit. However, with the re- 
alization of atomic gases with strong magnetic dipole mo- 
ments and the prospect of working with molecules that 
exhibit electric dipole moments [3, Q] , there is now sub- 
stantial interest in examining the physics of long-ranged 
interactions in these systems [3l ■ 

One important feature of long ranged interactions is 
that they can produce exceptionally intricate ground- 
state structures in the classical, or strong coupling, limit. 
A particularly elegant example of this is the case of clas- 
sical particles in a one-dimensional (ID) lattice interact- 
ing via an infinite-ranged convex potential studied by 
Pokrovsky and Uimin, and Hubbard (PUH) 0, i, Q. 
Here, it can be shown that the ground state filling frac- 
tion as a function of chemical potential /.t is a complete 
devil's staircase [loj . in which every rational filling frac- 
tion between and 1 is stable over a finite interval in /i, 
and the total measure of all such intervals exhausts the 
full range of ^. 

In this work, we show that the devil's staircase of 
PUH has dramatic consequences for the physics of quasi 
ID cold atomic gases. Building on the existing under- 
standing of this classical limit, we consider two pertur- 
bations of the devil's staircase that arise naturally in 
the experimental setting of cold atomic gases. The first 
of these is the introduction of a quantum kinetic en- 
ergy which now renders the problem sensitive to par- 
ticle statistics — which we take to be bosonic, given that 



this case can currently be realized with either atoms or 
molecules. The second perturbation involves tuning the 
onsite interaction independently of the rest of the in- 
teraction. This allows for a controlled departure from 
convexity, and hence from the PUH states considered 
previously. The first perturbation has the well under- 
stood effect of initiating a competition between the crys- 
talline, Mott phase that exists at zero hopping and the 
superfluid (Luttinger liquid in d = 1) that must exist 
at all fillings at sufficiently large hopping. By means 
of strong coupling perturbation theory similar to that 
previously studied in the Bose-Hubbard model [lH , and 
in extended Bose-Hubbard models with nearest-neighbor 
interactions l3, we derive a phase diagram that ex- 
hibits this evolution and which we supplement by stan- 
dard wisdom from the Luttinger liquid description of the 
transitions. The second perturbation introduces doubly 
occupied sites in the classical limit. While describing the 
resulting phase diagram in complete and rigorous detail 
is beyond the approaches we take in this paper, we give 
an account of the "staircase" structure of the initial in- 
stability and of regions of the phase diagram where the 
classical limit states exhibit superlattices of added charge 
built on underlying PUH states. At least some of these 
regions exhibit devil's staircases of their own. Finally, 
upon introducing hopping we are led to an infinite set of 
"supersolids" — which in this context are phases that are 
both Luttinger liquids and break discrete translational 
symmetries. 

Model: The Hamiltonian we will study is: 
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This model describes bosons in a deep ID optical lattice 
with hopping amplitude t, on-site interaction energy U , 
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and an infinite-ranged dipolc-dipole interaction Vo/r^. In 
cold-atom systems the hopping potential t is controlled 
by the depth of the optical lattice, and hence can be 
tuned over a wide range of values. In the cases where the 
bosons are atoms, e.g. ^^Cr the on-site interaction U 
may also be easily tuned using Feshbach resonances. We 
set the dipolar interaction to be maximally repulsive by 
aligning the dipoles perpendicular to the ID chain. 

If U is sufficiently large, the potential is everywhere 
convex [l^ and the t = (classical) ground states of ([T]) 
are those of a PUH Hamiltonian. For every rational fill- 
ing fraction u = p/q, the ground state is periodic with 
period q Q. Each such grou nd state is unique up to 
global lattice translations [14|. We denote these states 
commensurate ground states (CGS). Adding or removing 
a single particle from a CGS in the infinite volume limit 
produces a q-soliton state (qSS) containing q fractionally 
charged solitons of charge 1/q. Each soliton is a distor- 
tion of the periodic ground state which alters the distance 
between one pair of adjacent particles by 1. (For every 
commensurate state, there is a unique distortion of this 
type which minimizes the potential energy). The range 
of n over which each CGS is stable is given by 
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and falls off sharply as a function of g. At t = these 
intervals cover the entire range of fi pertinent to fillings 
less than unity, giving the devil's staircase structure. 

We first direct our attention to the t phase portrait 
but still at large U, where the potential is everywhere 
convex. The qualitative behavior of our system in this 
regime is reminiscent of the Bose-Hubbard model, with 
commensurate Mott lobes ceding to superfluid states as t 
increases. A state with large q may be treated as a state 
with smaller g at a nearby filling in which a crystal of di- 
lute solitons has formed. Hopping tends to liquify dilute 
crystals of solitons: at large separation the inter-soliton 
repulsion is smaller than the kinetic energy gained from 
delocalization. The delocalized solitons destroy long- 
ranged spatial order, creating a Luttinger liquid with full 
translational symmetry. Hence, as t increases, the system 
undergoes a transition from the Mott insulating CGS to 
a Luttinger liquid state, with larger q states liquifying at 
smaller t. 

Strong Coupling Expansion: To find the position of 
the phase boundary, we generalize the method of Ref . [Hi 
and compare the energies of the CGS and its adjacent 
qSS to third order in t using standard time-independent 
perturbation theory. This approach assumes that the 
phase transition is continuous, so that for a given t > 
0, values of fi for which Eqss(fJ') = Eccsdj) constitute 
the phase boundary. This assumption is well-founded, 
since the soliton repulsion ensures that the energy cost of 
creating multiple solitons is larger that of a single soliton, 



thus favoring a second-order transition. 

In a finite system one must account for the repulsion 
between solitons; here we drop these terms and consider 
the infinite volume limit. In the limit that the solitons 
are sufficiently well separated that we may neglect their 
interactions, the qSS is highly degenerate and can be ex- 
pressed in terms of a band of solitonic momentum eigen- 
states. Here we consider only the bottom of the band, 
which lies at zero momentum. 

The first and third order corrections to the CGS energy 
are zero. At v = p/q, the second-order correction is given 

by 



where lS.Ei ~ Ef*^ — E^'' is the difference in potential 
energies between the ground and the excited state formed 
by hopping from the i*^ occupied site in the ground state 
configuration. As the ground state is periodic, it suffices 
to calculate these energies for the p distinct particles in 
the repeated pattern. 

For the qSS at v ~ p/q, wc consider the infinite vol- 
ume limit, in which the total energy correction is simply 
q times the energy correction for a state with a single 
soliton. The soliton hops by q sites when a single boson 
on one of its edges is hopped by one site. The resulting 
energy corrections are 
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Here k is the soliton momentum, a is the lattice con- 
stant, and q is the denominator of the commensurate 
filling fraction, which appears here because one hopping 
displaces the soliton by q lattice spacings. The subscripts 
on Ej,-^±i indicate the distance between particle i and the 
soliton, and the direction of the hopping relative to the 
soliton. The special distance ri describes an excited state 
in which an anti-soliton is sandwiched between 2 solitons, 
or vice versa. These states contribute extra terms to the 
energy corrections of the qSS because of the ambiguity 
as to which soliton is associated with the ground state 
qSS. 

Fig. [1] shows the results of the perturbative calculation 
for selected filling fractions. The t ~ axis corresponds 
to the classical limit, in which the CGS states comprise a 
complete devil's staircase: every value of fi corresponds 
to a rationally filled ground state, except for a set of mea- 
sure 0. The figure shows the resulting Mott lobes: inside 
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FIG. 1: Perturbative calculation of the Mott to SF phase 
boundary in the plane, shown here for U = 20. Here 

the strength of all couplings is measured relative to that of the 
dipolar interaction strength Vb. Each lobe encloses a Mott in- 
sulating region in which the filling is fixed; the region outside 
the lobes is a superfluid of solitons. The 1/3, 1/5, 1/6, 2/7, 2/7 
and 3/7 -filled lobes are shown here. Every commensurate 
state of the complete Devil's staircase has a Mott lobe, but 
the range of hoppings over which a state exists falls off sharply 
with its denominator. 



each lobe the CGS is stable and the system is in a Mott 
insulating state. The Mott gap vanishes on the boundary 
of the lobe; outside of this region solitons proliferate and 
the system is in a Luttinger liquid phase. 

For any < > 0, only a finite number of insulating states 
exist; the rest are liquid states with a superfluid of con- 
densed solitons. The function is no longer a devil's 
staircase, but rather a piecewise smooth function, with 
plateaux of constant density separated by liquid phases 
whose density varies continuously with [i. The size of the 
commensurate region decreases sharply with q: states of 
higher g have both smaller ranges of stability in the clas- 
sical limit, and larger energy corrections relative to the 
CGS, so that the volume of the corresponding Mott lobe 
scales approximately as l/q^ [20| . The total volume occu- 
pied by liquid states can be estimated from the first-order 
approximation to the Mott lobe boundaries; we find that 
for small but fixed t the volume of the liquid region scales 
as approximately 1?!^ . At sufficiently large t we expect all 
insulating states to be unstable, and the particle density 
to vary smoothly with /i. 

Thus far, our results only apply to a homogeneous 
system, but we can estimate the effects of a harmonic 
trapping potential (present in current cold-atom exper- 
iments) using the local density approximation - this 
assumes that the trapping potential is slowly varying 
enough that it can be simply incorporated into the chem- 
ical potential, resulting in a spatially-dependent /x. Tra- 
jectories along the ID chain then correspond to cuts in 
Fig. [1] at fixed t/Vo. Thus, at t = we expect differ- 
ent commensurate fillings at different points along the 
trapped chain, with the most stable states (fillings 1/2 



and 1/3) occupying the largest regions within the trap. 
Note, however, that commensurate states with a period 
q greater than the lengthscale over which the trapping 
potential varies do not exist in the trapped system since 
they violate the local density approximation. Moreover, 
once t is nonzero but small, we have small islands of Mott 
states at various fillings p/g, separated by regions of su- 
perfluid which interpolate continuously between the two 
commensurate densities. 

The qualitative nature of the Mott transition can be 
deduced from existing knowledge of ID commensurate- 
incommensurate phase transitions, which we summarize 
here. Bosonization can be used to treat the kinetic term 
and dipolar interactions exactly; the lattice potential 
must be treated perturbatively. In bosonized form, our 
system is described by a Hamiltonian that is a sum of 
a Luttinger piece, accounting for the kinetic term and 
dipolar interactions [2l] |. and a sine-Gordon term which 
emulates the lattice in the continuum limit. 

The resulting Hamiltonian is well studied in both the 
context of the Mott transition [3l and the Frenkel- 
Kontorowa model of surface interfaces [l^. The upper 
and lower Mott lobes join in a cusp which is not accu- 
rately described by the perturbation theory- since 'small' 
hopping implies that the ratio t/5Ei of the hopping rel- 
ative to the energy gap to the solitonic states must be 
small, the range of t over which the perturbative treat- 
ment is valid decreases with q and never encompasses this 
point. The bosonized treatment reveals that crossing the 
edge of the Mott-Hubbard lobe induces one of two dif- 
ferent types of phase transitions. At the cusp joining the 
upper and lower Mott lobes, a constant-density phase 
transition of the Kosterlitz-Thouless type occurs. Every- 
where else the transition is well described by a simple 
two-band model with quasi-particles that are gapped in 
the commensurate phase, and with a density increasing 
as — /ic near the transition on the liquid side. 
Av^ray from the convex limit: The PUH CGS are the 
classical ground states so long as the potential is every- 
where convex. Since the on-site potential U is tunable 
experimentally, it is interesting to ask what happens to 
these states as U is lowered away from convexity and dou- 
ble occupancies begin to form. Of course, as U is lowered 
still further, triple and higher occupancies will also form 
but we will not push our analysis that far — the doubly 
occupied regime has enough challenges of its own. Our 
results on this are summarized in Fig. ^ and we now 
describe the analysis behind these. 

Let U[^^^\v) be the thresholds at which the PUH 
ground states give way to ones with at least one dou- 
ble occupancy — these are marked as the red points in 
Fig. (12). Observe that these thresholds increase mono- 
tonically with v, ui^^^\v') > ui^^^\v) for v' > v. 
This behavior can be qualitatively understood by con- 
sidering the implications of convexity at a given filling. 
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SufRcient conditions for convexity 0] are that, for all x, 

i(F(0) + Vi2x)) > Vix) and U > fvo (^^) (5) 

For convexity to hold everywhere, ^ must hold for 
a; = 1; below this threshold double occupancies may oc- 
cur. However, when perturbing about a given convex so- 
lution at fixed /x, solutions will be stable approximately 
until U violates ([5]) for x = rm, the minimal inter-particle 
distance. This implies that states with lower filling frac- 
tions are more stable against double occupancies, as the 
potential gain in lattice energy from doubly occupying a 
site is smaller. This is in contrast to the stability of the 
commensurate states as t increases, where the denomi- 
nator of the filling fraction determines stability. 

While our primary interest is in ground state transi- 
tions, much insight is gained by also computing a sec- 
ond threshold at a given filling. Past this threshold, 
U = Uc''^^\ a particle added to the state goes in as 
a double occupancy instead of fractionalizing into q soli- 
tons. From our above considerations, at fixed filling, the 
CGS is always more stable against forming double oc- 
cupancies than the qSS; a simple estimate suggests that 
the gap separating these instabilities is approximately 
the difference in energies for adding and removing a par- 
ticle to the state- which is equal to the interval in fi over 
which the state is stable at infinite U. 

Now consider a filling i/' > ly. As i/' can be constructed 
by adding charge to i', wc conclude that ui'~^'^^\v') > 
jj{qSS) ^^y^ the latter corresponds to the threshold at 
which adding an extra filling v' — v va. the form of the soli- 
tons of V loses out to adding it in the form of double oc- 
cupancies on top of the PUH state at v. As filling factors 
only slightly greater than v involve a dilute addition of 



charges, we further conclude that lim^,/ 



(CGS) 



jj(qSS) i^^y A somewhat more involved argument shows 

that lim,,^,- C/i^^^V) = Uf°^\v) H 

Our considerations so far imply intricate behavior for 
the location of the initial ground state instability, namely 
that Uc {v) is a monotone increasing function of v on 
the set of rationals with a discontinuity at each rational 
value of v. 

hm [/f > C/f ""^^M = lim 

We can understand the coarse features of this instability 
by considering first the most stable states (of denomina- 
tor q < qo, for some qo small enough to allow the exact 
thresholds to be computed), and deducing the expected 
behavior at fillings close to these. Fig. [2] plots Uc'^^^'^ 
(blue) and ui'~^'^^^ (red) for states with g < 15 in the 
vicinity of half-filling. The values of ui'^^^^ shown there 
are obtained by numerical minimization in the sector 
with one added charge at the specified filling. The values 



{u') 



(CGS) 

of ur ' are obtained by numerical minimization over 
configurations at the specified filling that contain exactly 
one double occupancy. This is the correct answer for all 
V = p/q for which p and q are not both odd. In such 
cases the double occupancy and its surrounding charge 
rearrangements give rise to even moments starting with 
a quadrupole. Consequently, double occupancies repel at 
all distances and enter via a continuous transition at the 
computed threshold. However, when p and q are both 
odd, the double occupancy has a dipole moment and the 
transition goes first order which causes the true C/i^^^^ 
to lie above our numerically determined value. This gap 
is not very large though, and vanishes at large q at least 
as 1/q'^ and we ignore it here. 
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FIG. 2: (color online) Numerically calculated values of 
jj{CGS) j^j.gj^^ gj^^ JJc'^^^'^ (blue), for a selection of filling frac- 
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tions u. Quoted values of U are measured relative to Vq. A 
segment of the j/-axis between Uc^'^^\l/2) and ui''^^\l/2) 
has been removed for better resolution of the rest of the phase 
diagram. The green and yellow boxes, bordered by black dots, 
indicate the approximate regions where the ground states con- 
sist of double occupancies in the 1/3 and 1/2-filled states 
respectively. Between the shaded regions the ground states 
consist of double occupancies on other, higher-denominator 
states. 

At any given filling, tracking the evolution of the 
ground state with decreasing U after double occupancies 
have been introduced is a problem of considerable com- 
plexity. Here we use the ideas developed thus far to iden- 
tify a family of regions in the (i/, U) plane where simpler 
descriptions emerge — these are indicated, in two simple 
cases, by the shaded regions on the figure. The basic idea 
is that once the v = p/q qSS becomes unstable to double 
occupancy, any particles added to the v = pj q state will 
be added as double occupancies, since these repel less 
strongly than solitons. Hence at first sight we expect, 
in the region ui'^^^\p/q) > U > ul^'^^\p/ q), states of 
filling V > p/q to consist of double occupancies in the 
v = p/q state. At rational fillings the double occupan- 
cies will arrange themselves in a crystal thus generating 
a commcnsuration distinct from that of the underlying 
p/q state — wc will refer to these as doubly-commensurate 
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states. However, for any given v this argument becomes 
less reliable as the density of added charge increases, since 
the background parent state itself becomes less stable to 
forming extra double occupancies. This can lead to tran- 
sitions in which the structure of the CGS collapses to a 
crystal of double occupancies over a background of sig- 
nificantly smaller filling. We have carried out a simple 
analysis of the location of this instability for the 1 /2 and 
1/3 plus double occupancy regions in Fig. ([2]) at selected 
fillings and these are marked by the black dots in the 
figure [23|. While this indicates that there are sizeable 
regions which can be described as simple descendants of 
the 1/2 and 1/3 states, we are not at present able to 
estimate the sizes of analogous regions for higher denom- 
inator fractions. Of course, to have a full solution of this 
would be equivalent to tracking the evolution of each v 
as U is decreased from Uc^'^^K 

As the barrier to triple occupancies is 3U , instabilities 
towards triple occupancy at a given filling will set in at 
approximately one third the value of U for instabilities 
to double occupancy. For U close to Uf ^{1/2), triple 
occupancies will not be favored at any filling fraction. 
A new staircase: The above discussion has led us to 
the doubly-commensurate states in the {f, U) phase di- 
agram: we remind the reader that such states are con- 
structed by periodically doubly occupying some fraction 
of the sites in a CGS. We now observe that in at least 
some cases these doubly commensurate states can form 
a devil's staircase of their own. 

First consider states constructed from double occu- 
pancies in the 1/2-filled state, which exist in the region 
shaded in yellow in Fig. ([2]). The energetics of such states 
can be divided into a) the constant interaction of the 
parent 1/2-fillcd PUH configuration with itself, b) the 
constant interaction of the added charges, irrespective of 
their location, with the parent 1/2 filled configuration 
and c) the interaction of the added charges with them- 
selves. This last part involves an interaction between the 
added charges which is convex again and thus leads to 
PUH configurations sitting on a lattice with a doubled 
lattice constant. The energy cost of adding a single dou- 
ble occupancy is U + Vd, Hence the doubly-occupied sites 
comprise a Devil's staircase with /i — > /j + U + Vd, and the 
widths of all intervals decreased by a factor of 8. Here 



is the interaction energy of each double 



occupancy with the underlying 1/2-filled state. At fixed 
U, this staircase is complete over the range of fillings 
for which increasing the particle density infinitesimally 
does not induce 'excess' double occupancies to form in 
the half- filled background lattice. In the case of the 1/2- 
filled state, for U sufficiently close to the upper cutoff 
this gives a complete staircase on 1/2 < < 1. Similar 
structures exist for all 1 / g- filled states in the appropriate 
range of U. As mentioned before, we do not, at present, 
understand the situation for doubly commensurate de- 
scendants of general rational fillings. 



"Supersolids": Thus far our considerations away from 
the convex limit have been purely classical. But we can 
equally consider states obtained from these modified clas- 
sical states upon the introduction of hopping. Specifi- 
cally, let us consider the fate of the doubly commensurate 
descendants of the PUH 1/q states considered above. 

In a manner entirely analogous to the problem with 
which we began this paper, the superlattice of added 
charges can melt via the motion of its solitons as t is in- 
creased resulting in a phase transition between the dou- 
bly commensurate state and a "super-solid" like phase 
in which the background 1/g-filled CGS coexists with a 
Luttingcr liquid. This is, in a sense the d = 1 version 
of the supcrsolid in higher dimensions, but it is worth 
noting that the d = 1 version in our problem exhibits a 
more divergent CDW susceptibility than superfluid sus- 
ceptibility as T ^ 0. 

To get a more quantitative account of these new 
phases, we may repeat the strong coupling treatment 
above. Fig. Q shows the phase portrait at intermedi- 
ate values of U near ly = 1/2. The black line traces 
the infinite-?/ Mott lobe, over which the background 
1/2-fillcd state is stable against forming solitons. The 
red line shows the threshold at which it is energeti- 
cally favorable to add a single double occupancy to the 
1/2 filled state. The blue curves show the positions 
of the Mott lobes for the doubly commensurate states. 
The presence of double occupancies stabilizes the 1/2- 
filled state against proliferation of solitons, so that the 
background remains commensurate at least within the 
infinite-C/ 1/2-filled Mott lobe, shown in black. This 1/2- 
filled super-solid phase has also been shown to exist in 
an extended Bose-Hubbard model with second-neighbor 
repulsion 17[ ; these numerical results are consistent with 
the phase portrait shown here for small t/V. 
Concluding Remarks: The long range of the dipolar 
interaction does produce, as promised, intricate phase 
diagrams for the one dimensional dipolar bosonic gas. 
Particularly striking are the singular staircase functions 
that showed up in our analysis in three different settings: 
in the i^ifi) curve for the exactly dipolar classical prob- 
lem, in the function uj^^'^^^ which marks the instability 
of the PUH states when the onsite U is tuned down and 
in the curves in selected regions of the {v^ U) plane. 
The other main set of results pertain to the presence of 
a large, indeed, infinite number of transitions between 
Mott crystals and Luttinger liquids or supersolids. The 
challenge of observing some of this physics in cold atomic 
gases is not trivial — the major obstacles are getting a rea- 
sonable simulacrum of a one dimensional gas of infinite 
extent. On the positive side, the control parameters we 
study here are eminently tunable and when combined 
with the physics we describe above, we hope, will moti- 
vate our experimental colleagues to invest some time in 
performing the needed miracles. 
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FIG. 3: (color online) Phase portrait in the vicinity of doubly 
commensurate states about 1/2 filling at f/ = Vo- Again, ^,t, 
and U in the figure are measured relative to Vo- The black 
lines indicate the boundaries of the Mott lobe ai U — 20. The 
red curve shows the chemical potential at which it becomes 
energetically favorable to add particles to the half-filled state 
as double occupancies, for U = Vo- The blue curves show 
some Mott lobes of the doubly-occupied staircase region; the 
region between these and the black lines (which delineate the 
region of stability of the 1/2-filled state at infinite U) is a 
super-solid state. 
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